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Abstract 

We studied possible evolution of the spin period and the magnetic field of the X- 
ray source RX J0720. 4-3125 assuming this source to be an isolated neutron star ac- 
creting interstellar medium. Magnetic field of the source is estimated to be 10^ — 10^ 
G, and it is difficult to explain observed spin period 8.38 s without invoking hypoth- 
esis of the magnetic field decay. We used the model of ohmic decay of the crustal 
magnetic field. The estimates of accretion rate (10~^^ — lO~^^M0/yr), velocity of the 
source relative to interstellar medium (10 — 50 km/s), neutron star age (2-10^ — lO^'' 
yrs) are obtained. 
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1 Introduction 

Isolated neutron stars (INS) have received special attention in the last few years. The idea 
of observing such objects in the X-ray range has emerged rather long ago (Ostriker et al., 
1970). Treves and Colpi (1991) supposed that INSs accreting matter of the interstellar 
medium (ISM) can be observed with the ROSAT satellite in the UV and X-ray ranges. 
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The estimations of energy characteristics, spectra and possibihties of observing INSs 
were made in several works (see, for example, Bohringer et al., 1987; Treves and Colpi, 
1991; Blaes and Madau, 1993). A large number of works is devoted to the study of the 
spatial distribution of INSs (see, for example, Blaes and Rajagopal, 1991). 

In the mid-1996, Haberl et.al. reported the discovery of the pulsating source RX 
J0720. 4-3125 with the ROSAT satellite in the soft X-ray range. We shall use two obser- 
vational characteristics of this source (Haberl et.al., 1996): the period p — 8.38 s and the 
blackbody temperature T = (79 ± 4) eV. 

Using the hypothesis that this source is an INS accreting matter from the ISM (we 
do not follow alternative hypotheses: further observations will show if it is true or not), 
we estimate the accretion rate and the magnetic field strength of this source. We show 
that the INS could not increase its period to the observed value over the Universe hfetime 
if we assume that it was born with the present-day magnetic field strength. Finally we 
conclude that the magnetic field of this neutron star (NS) had to decay. Then, using 
the model of ohmic dissipation of the magnetic field we calculate the magnetorotational 
evolution of the NS. 

The paper is organized as follows. In Section 2, wc present the reference material on 
the evolution of the neutron star period and on the model of ohmic dissipation of the 
magnetic field; in Section 3, the analytical estimates of the neutron star parameters are 
made; the evolutionary tracks of the INS are computed in Section 4; and the conclusions 
are given in Section 5. 

2 Evolution of the spin period and the magnetic field 
of the neutron star 

2.1 Spin period 

In the low-density plasma an INS may have four possible evolutionary states (Lipunov, 
1987): ejector, propeller, accretor, and georotator. A particular state is determined by 
the relation between four characteristic radii: the light cylinder radius Ri = c/uj; the 
stopping radius Rgt', the radius of the gravitational capture Rq = {2GM)/v'^; and the 
corotation radius Rco = (GM/cu^)^/^. Here, M is the NS mass, c is the speed of hght, to 
is the spin frequency, and v^o is the NS velocity with respect to the ISM. 

The relationship between these radii determines two critical periods: Pe and Pa which 
separate different stages of the NS evolution. These periods can be estimated using the 
formulae (Lipunov, 1987): 
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Here /i is the magnetic dipole moment, M = ttRqPVoo is the accretion rate, p is the ISM 
density, and kt is a dimensionless constant of the order of unity. 

If p < Pe, then NS is at the ejector stage; if < p < P^, we have the NS at the 
propeller stage; and if p > Pa and Rgt < Rg: then NS is an accretor. In some cases the 
situation is possible when p > Pa, but Rst > Rg and accretion is impossible because of 
the formation of the geolike magnetosphere. However, we shall not be interested in the 
georotator stage since we consider the accreting NSs for which Rst < Rg- 

At the ejector stage, the evolution of the spin period is determined by the losses of 
the INS kinetic energy due to magnetic dipole radiation: 

P^'-^-^, (3) 

where R is the NS radius, / is the moment of inertia, and B — p./ R? is the magnetic field 
strength. 

At the propeller stage, the NS spin-down rate is determined by the transfer of the 
angular momentum to the surrounding matter (lUarionov and Syunyaev, 1975): 



p = p^B^"M'l\ (4) 

TT I 
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At the accretion stage, the NS is acted upon by two moments of forces: 

d(27rl/p) 



dt 



Ksd + Kturb, (5) 



R^ 

Here Kgd is the braking moment of forces which can arise due to the possible turbulization 
of the ISM. Kturh acts randomly and can either spin up or spin down the NS (Lipunov 

and Popov, 1995). 

The change in the period of an accreting INS is due to its interaction with the turbu- 
lized ISM. This introduces specific features into the problem of the period evolution. If 
we adopt the hypothesis of spin acceleration of the NS in the turbulized ISM (Lipunov 
and Popov, 1995), the new characteristic period emerges, 

Pe, = 960^yV^^^4/V:^f -S^-^''^^4'' sec = (6) 
= U^^kr^,%'ll^MZl^vyiv;,'''M;r sec. 
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where /X3o is the magnetic dipole moment in units of 10^° Gs cm^, I45 is the moment of 
inertia in units of 10^^ g cm^, p_24 is the ISM density in units of 10~^^ g cm~^, Vooq is the 
NS velocity relative to the ISM in units 10^ cm/s, and Vt^ is the turbulent velocity in units 
of 10^ cm/s. The period Pgg corresponds to the NS rms rotation rate obtained from the 
solution of the corresponding Fokker-Planck equation. In reahty, the rotational period of 
INS fluctuates around this value. Note that here we make a more accurate estimate of 
the period than that made in the work of Lipunov and Popov (1995) (we are grateful to 
M.E. Prokhorov for his assistance in performing the corresponding calculations). We take 
into account the three-dimensional character of turbulence, i.e. the fact that the vortex 
can be oriented not only in the equatorial plane but also at any angle to this plane. In 
this case, diffusion occurs in the three-dimensional space of angular velocities. 

2.2 Ohmic dissipation of the crustal magnetic field of the neu- 
tron star 

The magnetic field decay in the neutron star crust has been investigated by many authors 
(see, f.e., Urpin and Muslimov, 1992). Such a decay is governed by the induction equation 

where a is the conductivity; v is the velocity of the crust motion; and v = at the ejector 
and propeller stages. At the accretor stage in spherically symmetrical case, v = {—Vr, 0, 0), 
with the component 

M 

where r is the distance to the NS center, p{z) is the density of matter at the depth 
z — R — r. Equation (7) is highly simphfied in the case of the dipole field and reduced to 
the one- dimensional parabolic equation for the vector potential. The boundary conditions 
are set in the same ways as in the work of Urpin and Muslimov (1992). 

The conductivity a is determined basically by the scattering of electrons by phonons 
and impurities: 

11^1 

The phonon conductivity cTph depending on density and temperature dominates at 
high temperatures and at moderately high densities. At lower temperatures and higher 
densities the impurity conductivity (Jimp prevails. The impurity conductivity is indepen- 
dent of the temperature but is the function of concentration and the charge of impurities 
which are characterized by the parameter 
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where n and Z are, respectively, the concentration and the charge of the major ion species, 
and n' and Z' are the concentration and the charge of the impurity, the summation is over 
all impurity species. The analytical formula for the impurity conductivity was obtained 
by Yakovlev and Urpin (1980), while the phonon conductivity was taken from the work 
of Itoh et al. (1993). The initial field is assumed to be concentrated mainly in the surface 
layer of a certain thickness. The density of matter, po, at the inner boundary of this layer 
is the parameter of the model. The quantity Q is assumed to be independent of the depth. 

The main results of the calculation of the decay of the INS magnetic field are the 
following. The magnetic field dissipation proves to be intimately related to the NS thermal 
evolution. For the standard cooling, at which the NS neutrino luminosity is determined 
chiefly by the modified URCA-processes (Pethick, 1992), in the first million years the field 
decays by a factor of 2 — 1000, depending on the initial depth of the current layer and on 
the adopted equation of state for the star core (Urpin and Koncnkov, 1997). As the NS 
cools down, the conductivity increases and the field decay slows down. The decay rate 
at the later stage depends on aimp and, therefore, on Q. For example, when Q — 0.01, 
the field does not decrease in the subsequent 10^ yrs. However, as soon as the magnetic 
field diffuses through the entire crust and reaches the superfluid core (after ~ 10^ yrs at 
Q = 0.01). the decay becomes exponential. 

Figure 1 shows the decrease of the NS surface magnetic field with time for various 
parameters po and Q. We perfomcd our calculations for the model of the NS based on 
the Friedman-Pandharipande (1981) moderately hard equation of state in the star core, 
with the mass of the neutron star M = IAMq, the radius R = 10.6 km, and the crust 
thickness AR — 940 m (Van Riper, 1988). It is apparent that the decay at the initial 
stage {t < 10^ yr) is determined by the initial depth of the current layer and that the 
subsequent decay rate depends on the amount of impurities. 

Accretion may affect the field evolution. First, it heats the NS crust (Zdunik et al., 
1992) and, therefore, decreases the conductivity. Second, the fiux of matter toward the 
star center arises; it tends to bring the field to deeper layers. However, calculations 
show (Urpin et al., 1996) that accretion with the rate M < lO~^"^M0yr~^ speeds up 
insignificantly the field decay. 

3 Estimates of parameters of RX J0720.4-3125 

Using the hypothesis that we observe an accreting isolated NS and taking the observed 
period p — 8.38 s and the temperature T — (79 ± 4) eV (Haberl et. al., 1996), we can 
obtain the constraints on the magnetic field, accretion rate, and luminosity. 
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Figure 1: The change of the surface magnetic field of the isolated neutron star with time 
for the model of standard cooling. Curves 1, 2, and 3 correspond to the initial depths of 
the current layer, 10^^, 10^^, and 10^'^ g cm~^, respectively. The solid curves correspond 
to Q = 0.001; the dashed curves, to Q = 0.01; the dot-dashed curves, to Q = 0.1. 
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1. If we really observe the accretor, then p > Pa{B, M), where Pa is defined by (2). 
This yields the first constraint on B and M: 

S < 4.6 ■ IO^M^^^R-^^M]/^^ G, (8) 

2. We see the pulsating radiation, it means that the accreting matter is funneled by 



the magnetic field onto the polar caps. Using the known relation i?cap = \J{R/Ra) ■ R, 
connecting the polar cap radius Reap, the Alfven radius Ra, and the NS radius, we can 
estimate Reap as 

i?eap = 0.53 • B^'/'m'SIrT'mI!:' km, 

where Bg — G. Prom the condition Reap < R (or Ra > R) we obtain another 

constraint: 

S > 3.4 • 10^Mi//i?6 '/^M^is G. (8) 
3. And finally, we know the polar cap temperature. Prom the condition 

^ „2 ^4 GMM 

cap ^ 

we find the additional relation between B and M: 

5 = 2.5- 10^ (r/79 eV)^ RIM^^^'"^ MZlt G. (9) 

Combining inequalities (8) and (9) and condition (10) yields the allowed range of the B 
and M values: 

2 • 10^ < B[G\ < 10^ 
6 ■ 10-1^ < M[Mo/yr] < 4 ■ IQ-^l (10) 

This range of possible M values corresponds to the allowed range of luminosities L: 

7 ■ 10^^ < L[erg/s] < 5 ■ 10^^ 

Let us estimate the spin-down time to p = 8.38 s assuming that the NS was born as a 
radio pulsar with po <^pe- The time of spin-down to p = p£; at a constant magnetic field 
is deduced from (3) and (1): 

_ 3c3/ „ _ 



4 • lQ''MZll''h,v^i'RfB^' yrs 
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For B < 10^*^ G and M < 10~^^MQyr~^ the spin-down time exceeds the age of the 
Universe. The spin-down at the propeller stage takes an additional time. Thus, we 
come to the conclusion that with the estimated values of the magnetic field strength, the 
accretion rate, and the velocity of motion relative to the ISM the NS could not spin down 
to p = 8.38 s. This means that in the past the given INS had a larger magnetic field that 
has decayed over the time of evolution. 

4 The evolutionary tracks of neutron star on the B — 
P diagram 

The evolution of the spin period was calculated from formulae (3) and (4) taking into 
account the decay of the magnetic field at the ejector and propeller stages. Since the 
treatment of the spin-up rate in the turbulized ISM at the accretor stage is not quite a 
simple matter, we applied to the description of the period evolution the following simpli- 
fying model. When the NS gets to the accretor stage, the spin-down momentum substan- 
tially exceeds the spin-up momentum because the acceleration occurs in the turbulized 
ISM and there is no constant spin-up momentum, unlike in a binary system. However, 
we can obtain an analog of the equilibrium period (Section 2.1) which corresponds to the 
stationary solution of the Fokker-Planck equation (Lipunov, 1987; Lipunov and Popov, 
1995). Therefore the evolution of the rotational period at this stage was treated as a 
stationary spindown from to Peg after which the period was set to be Peg, which, in 
its turn, changed due to the magnetic field decay. 

We calculated the magnetic and spin evolution of the NS with mass M = IAMq 
for the accretion rates lO~^^M0yr~^ and 10~^^ MqJT^^, moving through the ISM with 
density p = lO^^'' g/cm^. Such accretion rates at the given ISM density correspond to 
the velocities of motion 19 and 41 km/s. We assumed the NS to be born as an ordinary 
pulsar with a short period (po — 0.01 s) and a "standard" initial magnetic field Bq — 10^^ 
G. The model of the ohmic decay of magnetic field allows one to obtain both the high and 
low dissipation rate. However, in this case we should fit the parameters of the model so as 
to obtain the accreting INS with a period of 8.38 s and a field of 2 • 10^ G for the accretion 
rate M = IO^^^Mq yr~^ or with a field of 4 • 10^ G for the accretion rate lO^^^M© yr~^ 
(10). Therefore, the field should decay with a moderate rate to give the NS enough time 
to slow down to p = 8.38 s. Moreover, it is desirable at the initial stage of evolution 
to provide the agreement with magnetic fields and periods of the observed radio pulsars. 
Basing on these assumptions we can obtain the estimates of the parameters po and Q, as 
well as the lower hmit on the INS age. 

The evolutionary tracks of the NS for the accretion rates M — IO'^^Mq yr~^ and 
M — IO^^^Mq yr'^ are shown in panels a) and b) of the Fig. 2. 
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Figure 2: The evolutionary tracks of the neutron star for the accretion rates M = 
lO^^^M0yr~^ (a) and M = IO^^^Mq yr~^ (b). The model parameters are described in 
the text. The dashed hnes correspond to p = P^; the dot-dashed lines, to p = Pa- The 
dashed line in Fig. 2a shows for the second track the neutron star evolution with no 
acceleration in the turbulized intestellar medium. The numbers near the marks in tracks 
denote the logarithm of the neutron star age in years. The observed radio pulsars are 
indicated by dots (Taylor et al., 1993). 
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Tracks 1 in both panels of Fig. 2 illustrate the evolution with the maximum possible 
rate of magnetic field dissipation with time. The model parameters for the accretion rate 
M = IO-^^Mq yr-i are: po = 3 ■ 10^^ g/cm^ Q = 0.02, and v = 10^ cm/s. The NS was 
born at the point (po, Bq) at the left upper corner of the B — P diagram. In the first lOg 
yrs, the field decays by about 20 times, while the NS slows down top 0.3 s. At that time 
the NS represents a typical radio pulsar. Subsequently, due to the cooling, the field decay 
slows down. In the next 8 ■ 10'^ yrs the field decays to 4.5 ■ lO"*^^ G, the period increases, 
and the radiopulsar dies. However, the ejector stage continues for 1.2 ■ 10^ years. During 
this stage, the field decays to 2.5 • 10^° G, the period increases to 2.1 s, and the NS goes 
over to the propeller stage (1) which lasts ~ 10^ yrs. At this stage, further deceleration 
of the NS occurs, according to (4), to p = 5.7 s, and the field decays to 3 ■ 10^ G. The 
star goes over to the accretor stage and becomes the source of periodic X-ray radiation. 
The core heating speeds up insignificantly the field decay (Urpin et.al., 1996). The field 
decays to 2 ■ 10^ G over 4 ■ 10 yrs. The period does not increase because the field is weak 
enough. As soon as Peq becomes equal to the current period, the spin-up of the NS due 
to the interaction with the ISM may occur. The period fiuctuates around Peq, but we do 
not estimate here the amplitude of these fluctuations. 

Tracks 2 illustrate the evolution with slower field decay. Such a decay can be obtained, 
for example, for the crust with lower impurity content. For M — 1Q~^^Mq yr~^, we set 
Q = 0.01, while the other parameters remained the same as for track 1 in Fig. 2a. Tracks 
1 and 2 coincide at the initial stage of evolution, when the dissipation rate does not depend 
on the impurity concentration. However, in 10^ yrs the tracks diverge. As a result, the 
star resides at the ejector stage 2 • 10^ yrs going over to the propeller stage with p = 3 
s and B = 4 ■ 10^° G. The spin-down rate at the propeller stage is much higher than in 
the first case for three reasons: due to the longer spin period, the stronger magnetic field 
during the NS transition from the ejector to the propeller stage, and to the lower rate of 
the field decay caused by the lower value of the Q parameter. At the accretor stage, the 
NS spin period decreases to 190 s over 2-10^ yrs, and when the field decays to 4-10^ G, the 
effect of acceleration in the turbulized ISM may begin to act. The track in the absence of 
such acceleration {vt — 0) is shown by the dashed line. In this case the final period turns 
out to be about 200 s. The turbulent acceleration may decrease this period. However, 
the self-consistent calculation of the evolution of the NS period at the accretor stage with 
allowance for the magnetic field decay would require the solution of the Fokker-Planck 
equation for the distribution function in the space of angular velocities, which is beyond 
the scope of this work. 

At the accretion rate M — 1O~^^M0 yr~^ we used for track 1 the following parameters: 
Po = 3 ■ 10^^ g/cm^, Q — 0.01, and Vt — 10^ cm/s. The NS comes to the propeller stage 
in 2.7 ■ 10^ yrs having a period of 3.1 s and a magnetic field of 2 ■ 10^° G. The transition 
to the accretor stage occurs when the NS has a period of 4.9 s and a magnetic field of 



10 



7 • 10® G. Deceleration at this stage does not occur because the magnetic field is low. The 
full time taken for the field decay to 4 ■ 10^ G is 6 • 10^ yrs. In this case p = 5 s. The 
second track in Fig. 2b was calculated for the greater initial depth of the current layer 
corresponding to po = 6 ■ 10^^ g/cm^. As a consequence, at the initial stage of evolution 
the magnetic field decays with a slower rate, the ejector stage lasts 2.1 • 10^ yrs, and the 
propeller stage lasts 1.9 ■ 10^ yrs. For the field 4 ■ 10^ G, the period is 63 s. 

5 Conclusion 

The observed period and the temperature of the X-ray source RX J0720.4- 3125 can be 
accounted for by using the hypothesis for the ISM accretion onto an old INS. We showed 
that the NS magnetic field is low in this case {B < 10^ G). The time taken for deceleration 
to p = 8.38 s with such a magnetic field exceeds the age of the Universe. We supposed 
that the NS was born with higher value of magnetic field and that the field strength has 
substantially decreased during the evolution. Using the model of ohmic dissipation of the 
magnetic field in the NS crust, we calculated the possible evolution of the NS on the B — P 
diagram. The observed rotational period can be obtained at Q = 0.01 — 0.05. However, 
the evolution of the period depends on the field dissipation rate and, therefore, on the 
parameters of the decay model. Thus, the twofold change in the impurity parameter Q 
caused the spin period to change more than by the order of magnitude at the accretor 
stage (Fig. 2a, tracks 1 and 2). For this reason, observations of periods of old INSs may 
become an important test for the validity of the model of evolution of neutron stars. 

The decay of the magnetic field may affect the estimate of the total number of the 
observed accreting INSs. In particular, because the formation of a periodic X-ray source 
at the propeller stage calls for rather strong magnetic field (Popov, 1994), its decay may 
reduce the number of sources of this kind. At the same time, the absence of pulsations 
for other INS candidates may suggest that their field has already decayed to such extent 
that it cannot funnel the plasma toward the INS polar caps. 
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